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[i] Phytoplankton absoiption spectra and High-Performance Liquid Chromatography 
(HPLC) pigment observations from the Eastern U.S. and global observations from 
NASA’s SeaBASS archive are used in a linear inverse calculation to extract pigment- 
specific absorption spectra. Using these pigment-specific absoiption spectra to reconstruct 
the phytoplankton absoiption spectra results in high correlations at all visible wavelengths 
(r 2 from 0.83 to 0.98), and linear regressions (slopes ranging from 0.8 to 1.1). Higher 
correlations (r 2 from 0.75 to 1.00) are obtained in the visible portion of the spectra when 
the total phytoplankton absorption spectra are ‘unpackaged ’ by multiplying the entire 
spectra by a factor that sets the total absorption at 675 nm to that expected from absorption 
spectra reconstruction using measured pigment concentrations and laboratory-derived 
pigment-specific absorption spectra. The derived pigment-specific absorption spectra were 
further used with the total phytoplankton absoiption spectra in a second linear inverse 
calculation to estimate the various phytoplankton HPLC pigments. A comparison between 
the estimated and measured pigment concentrations for the 1 8 pigment fields showed good 
correlations (r 2 > 0.5) for 7 pigments and very good correlations (r 2 > 0.7) for chlorophyll a 
and fiicoxanthin. Higher correlations result when the analysis is carried out at more local 
geographic scales. The ability to estimate phytoplankton pigments using pigment-specific 
absorption spectra is critical for using hyperspectral inverse models to retrieve 
phytoplankton pigment concentrations and other Inherent Optical Properties (IOPs) from 
passive remote sensing observations. 

Citation: Moisan, J. R., T. A. H. Moisan, and M. A. Linkswiler (2011), An inverse modeling approach to estimating phytoplankton 
pigment concentrations from phytoplankton absorption spectra, J. Geophys. Res., 116, C09018, doi:10.1029/2010JC006786. 


1. Introduction 

[ 2 ] The light absorption properties of marine phytoplank- 
ton influence the manner in which solar energy radiates 
through the ocean and control the level of energy made 
available to phytoplankton for primary production. Phyto- 
plankton can absorb light across the visible and into the UV 
portion of the light spectrum. Knowledge on the shape of 
phytoplankton absorption spectra is a requirement in present 
inverse models that estimate phytoplankton chlorophyll 
concentrations [ IOCCG , 2006] and for input into bio-optical 
models that predict carbon fixation rates for the global ocean 
[Behrenfeld and Falkowski, 1997; Carr et al., 2006], The 
shape and the magnitude of the phytoplankton absorption 
spectra is controlled primarily by the concentration of various 
photosynthetic and photoprotective pigments and by the level 
of pigment package effect within the cells [Stuart etal., 1998; 
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Lohrenz et al, 2003], though the specific influence of these 
two processes varies with depth, phytoplankton species 
composition, cell size and physiology. 

[ 3 ] Earlier attempts at reconstructing phytoplankton in 
vivo absorption spectra from pigment concentrations used in 
vitro pigment-specific absorption measurements of pure 
pigment standards ‘shifted’ spectrally to match observed 
shifts in the in vivo absorption maxima positions from pig- 
ment protein complexes within the cell [Bidigare et al., 
1990]. This technique assumes that the reconstructed 
absorption spectra are a linear combination of the pigment- 
specific absorption spectra from individual “unpackaged” 
pigments, such that the in vivo, “unpackaged” phytoplankton 
absorption spectra 

m 

4*(A) = ^C,flr(A), (1) 

1=1 

where af( A) is the individual weight-specific ‘shifted’ pig- 
ment absorption spectra derived from absoiption spectra of 
individual pigments within a known solvent and C, is the 
concentration for the ith pigment. The terms ‘weight- 
specific’ and ‘pigment-specific’ absorption [m 2 mg *] are 
used interchangeably in the scientific literature and hold the 
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Figure 1. Weight-specific (or pigment-specific) in vitro 
absoiption spectra of various pigments, a* (A), derived from 
measuring the absorption spectra of individual pigments in 
solvent and shifting the maxima of the spectra according to 
Bidigare et al. [1990], Data obtained courtesy of Annick 
Bricaud [see Bricaud et al., 2004], Note that only 14 of the 
total 18 pigments encountered in this study have pigment- 
specific relationships available from laboratory studies. 


same definition of absolute absorption [m 1 ] per concentra- 
tion of a given pigment [mg m 3 ], 

[ 4 ] Nelson et al. [1993], using a similar combination of 
absorption spectra and HPLC pigment measurements, noted 
that when pigment package effects were small the Bidigare 
et al. [1990] reconstruction method compared well to the 
measured absorption and that underestimates from the 
reconstruction method was likely due to absorption from 
phycobiliproteins or other pigments not measured in tradi- 
tional HPLC analysis. 

[5] In addition to the indirect reconstruction method 
(equation (1)), Hoepffner and Sathyendranath [1991] 
developed a technique that fitted the phytoplankton 
absorption spectra to combinations of Gaussian curves/ 
spectra such that the reconstruction of the in vivo phyto- 
plankton absorption spectra 


L 

a pi: (A) = A ^ c/(3 maX) / exp 

i= 1 


(A - A,-) 2 ' 
2aj 


( 2 ) 


where L denotes the number of total Gaussian curves to be 
summed over, c, is the concentration of the pigment asso- 
ciated with the individual Gaussian curve, Umax,; is the 
pigment-specific maximum absorption level associated with 
the center of the Gaussian curve located at the A, : wave- 
length, and a i denotes the half width of the Gaussian curve. 
Analysis of the results from fitting 1 1 such Gaussian curves 
to a variety of phytoplankton culture absorption spectra from 


9 species within 3 major algal groups (3 Bacillariophyta, 
4 Chlorophyceae, and 2 Prymnesiophyceae) showed that 
Gaussian curve summations reconstruct the absorption 
spectra very well. However, attempts to model the vari- 
ability in the various Gaussian curve parameters as a linear 
function of specific pigment concentrations showed mixed 
results, with coefficient’s of determination, r 2 , ranging from 
0.33 to 0.99. Lutz et al. [1996] demonstrated that these 
Gaussian reconstructions of absorption spectra are improved 
if the ocean regions are partitioned according to their 
absorption characteristics. Using this technique, analysis of 
Arabian Sea and Vancouver Island seawater samples by 
Stuart et al. [1998] showed that 29^-2% of the variability 
in absorption at 440 nm was due to changes in pigment 
composition and that 58-71% of the variability was due to 
pigment package effects. In a similar study using seawater 
samples from the coasts of North Carolina and Florida, 
Lohrenz et al. [2003] also noted that pigment package 
effects accounted for up to 62% of the variability in the 
absorption spectra maximum amplitude bands (a* ax ,,j and 
that variations in pigment composition only accounted for 
10-28% of variations. 

[6] Further attempts at refining this method [Evans and 
Cornford, 2003] have expanded the linear fit approach for 
modeling the pigment-specific maximum absorption such 
that the reconstructed in vivo phytoplankton absorption 
spectra 


m L 

ftph (-^) = ^ ^ Q ^ ^ ^ma x,ij 
*=1 7=1 


(A -A,f 
2aj 


(3) 


such that a suite of L Gaussian curves (each with a unique 
maximum absorption a* aXjIi/ -, half width <jy and central 
wave band A /y ) are linked to n pigments. The utility of 
this approach is that the resulting equation set can be used 
to compute pigment-specific absorption spectra similar to 
those developed initially by Bidigare et al. [1990]. How- 
ever, implementation of this method has proved more 
challenging, likely due to correlations between pigment 
package effects and concentrations (D. Comford, personal 
communication, 2010). In these latter two methods, the 
derived pigment-specific absorption spectra are influenced 
by the level of pigment package effects that occurs within 
various samples. Also, the numerical retrievals are based 
upon individually solving the equations for each absorption 
sample, i.e., absorption values at different wavelengths from 
one sample are used to solve for the various unknowns in 
the equation. 

[ 7 ] A refinement of the Bidigare et al. [1990] recon- 
struction approach was carried out by Bricaud et al. [2004] 
through an expansion of the number of pigment-specific in 
vitro absorption spectra (Figure 1). Using these expanded 
number of spectra and HPLC pigment measurements and an 
independent assessment of the package effect index, the 
results showed that the dominant cause of deviations from 
the average absorption to chlorophyll a relationships is 
caused by cell size distributions in the phytoplankton pop- 
ulation, the primary factor controlling the pigment package 
effect. 

[8] Regardless of the reconstruction methods previously 
used, it remains difficult to separate out the influence of 
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Figure 2. Map of the (a) sample locations in the Maryland and Virginia, Gulf of Maine and Cape Cod 
coastal regions, and (b) sample locations from these samples combined with the NASA SeaBASS samples. 


pigment variability with pigment package effects since the 
pigment package effect is overall linked to cell size dis- 
tributions and therefore phytoplankton species composition. 
This link in spectral absorption and dominant cell size 
relationship has been well demonstrated by Ciotti et al. 
[2002] in an attempt to use the shape of the phytoplankton 
absorption spectra to estimate the dominant cell size. In 
addition, Moisan and Mitchell [1999] have shown that the 
pigment package effect for a single species can cause a 
greater than twofold variability in phytoplankton absorption 
spectra due to changes in light and temperature alone. 

[ 9 ] The need to understand the causes of variability in 
phytoplankton absorption spectra as it relates to phyto- 
plankton pigment composition, cell size distribution and 
pigment package effect is obvious. Present limitations in our 
knowledge in this arena have placed restrictions on usage of 
passive remote sensing reflectance spectra from satellite, 
aircraft or stationary platforms. Present-day ocean inverse 
models [IOCCG, 2006] utilize chlorophyll-specific absorp- 
tion spectra models, derived from either simple models such 
as a Gaussian model [Hoge and Lyon, 1996, 1999] or from 
inverse models \Maritorena et al., 2002], that are limited at 
being able to only retrieve a single phytoplankton pigment, 
chlorophyll a. The development of the next generation of 
ocean color satellites capable of hyperspectral observations 
(typically regarded as 5 nm or less spectral resolution 
between 400 nm to 700 nm) will allow us to make inroads 
into estimating additional phytoplankton pigments and 
possibly functional types and size spectra. 

[ 10 ] The objective of this paper is to analyze a novel and 
direct method of using multiple measurements of phyto- 
plankton absorption spectra and pigments to derive pig- 
ment-specific absorption spectra that can be used to directly 


estimate phytoplankton pigment concentrations from phy- 
toplankton absorption spectra. The work utilizes a range of 
absorption spectra and HPLC pigment measurements that 
have been collected in both Case I and Case II off the coasts 
of Maryland, Virginia, the Gulf of Maine, and Cape Cod 
(Figure 2a and Table 1.) and supplemented with additional 
observations obtained from the NASA SeaBASS data 
archive (Figure 2b). The methodology presented here can be 
incorporated directly into hyperspectral inverse models and 
employed for estimation of phytoplankton pigment con- 
centrations. Approaches such as these, which link phyto- 
plankton absorption spectra to pigment concentrations, are 
supportive of applications of future ocean color remote 
sensing satellite platforms (such as NASA’s Geostationary 
Coastal and Air Pollution Events [GeoCAPE], Aerosol- 
Clouds-Ecosystem [ACE], and Pre-ACE [PACE] missions) 


Table 1. Related Cruise and Sample Infomiation 


Cruise 

Location 

Dates 

Number of 
Samples 

15 COBY cross 
shelf surveys 

VA/MD cross 
shelf survey 

Monthly from 

101 

BIOME 3 

DE/MD/VA 
coastal region 

26-30 July 2005 

26 

BIOME 5 

DE/MD/VA 
coastal region 

9-12 May 2006 

19 

MAA 1 

Gulf of Maine 

26-30 April 2007 

38 

MAA 2 

Gulf of Maine 

26-28 May 2007 

50 

MAA 3 

Gulf of Maine/ 
Martha’s Vineyard 

6-8 June 2007 

32 

SeaBASS Archive 
Total Samples 

Global Extent 

Various 

608 

874 
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Table 2. List of Pigments Used in Analysis 


Number 

Name 3 

Field 

Names 3 

Pigment 

Classification 

1 

Chlorophyll c (cl + c2 + c3) 

chl c 

Photosynthetic 

2 

Chlorophyllide 

chlide 

Degradation Product 

3 

Pheophorbide 

phide 

Degradation Product 

4 

Peridinin 

peridinin 

Photosynthetic 

5 

19' B utanoyloxy fucoxanthin 

but-fuco; 
1 9 '-but 

Photosynthetic 

6 

Fucoxanthin 

fuco 

Photosynthetic 

7 

Neoxanthin 

neo 

Photosynthetic 

8 

Violaxanthin 

viola 

Photosynthetic 

9 

1 9' Hexanoyloxyfucoxanthin 

hex-fuco; 
19 '-hex 

Photosynthetic 

10 

Diadinoxanthin 

diadino 

Photoprotective 

11 

Alloxanthin 

alio 

Photosynthetic 

12 

Diatoxanthin 

diato 

Photoprotective 

13 

Zeaxanthin 

zea 

Photosynthetic 

14 

Lutein 

lut 

Photoprotective 

15 

Chlorophyll b 
(mono and di-vinyl) 

chl b 

Photosynthetic 

16 

Chlorophyll a 
(mono and di-vinyl) 

chl a 

Photosynthetic 

17 

Pheaophytin a 

phytin 

Degradation Product 

18 

Carotenoids 

car 

Photosynthetic and 
Photoprotective 


“Pigment names and field names obtained from SeaBASS Bio-Optical 
Archive. 


that will be hyperspectral in nature and allow for the full 
inversion of ocean color remote sensing reflectance. The 
added value of hyperspectral wavelength resolution is to 
allow for the comprehensive description of the phyto- 
plankton community pigment suite. These estimates could 
potentially be further exploited using applications such as 
CHEMTAX [Mackey et al., 1996] to estimate phytoplank- 
ton community composition. 

2. Methods 

2.1. Regional U.S. East Coast Absorption Spectra 
and HPLC Pigment Data Set 

[n] A total of 266 water samples were collected during 
20 different cruises from 2 U.S. eastern coastal ocean regions: 
1) Delaware, Maryland and Virginia (149 samples), and 2) the 
coastal waters within the Gulf of Maine and near Martha’s 
Vineyard (120 samples). The data set is inclusive of all 
seasons for the DL/MD/VA region, and of entirely spring 
samples from the Gulf of Maine (Table 1). All samples were 
gathered and processed utilizing similar methods. Water 
samples were collected from Niskin bottles at various depths, 
filtered through duplicate Whatman GFF filters, placed in 
Flistoprep capsules, and flash frozen in liquid nitrogen for 
later laboratory analysis. The majority of the samples were 
collected in the upper few meters of the water column, 
a region with the greatest influence due to phytoplankton 
pigmentation on passive radiance sensors, such as on-board 
aircraft or satellites. 

[ 12 ] Absorption Spectra. Phytoplankton absorption spec- 
tra samples were processed using the filter pad technique 
that partitions the particulate and detrital fraction [Kishino 
et ah, 1985] to yield a phytoplankton absorption coefficient 


[Mitchell, 1990]. Absorption spectra were acquired on a 
Perkin Elmer LS800 UV/VIS Spectrophotometer at 1 nm 
intervals from 300 nm to 800 nm using a 4 nm slit-width to 
fully resolve the spectral shape. The instrument was back- 
ground corrected before each batch of six samples, and 
baselines were run on filtered seawater. Initial samples were 
run using a filtered-seawater blank for comparison to deter- 
mine the absorption by all particles, a p (A). After particulate 
absorption measurements were carried out, the pigments were 
extracted by thrice placing the filters into 10 mL of 100% 
methanol and then rinsed with filtered seawater to remove 
phytoplankton pigments. The filters were then rescanned in 
the spectrophotometer, again using a filtered seawater blank 
for comparison, to determine the non-pigmented, or detrital 
absorption aj, A). Samples were extracted in methanol for 
a final 10-min period if peaks in absorption spectra appeared 
by the operator to be derived from residual pigments 
remaining in the sample. Several Beta-factors were compared 
using the data set from the NASA sponsored absorption 
workshop [Mitchell et ah, 2000], The Mitchell [1990] cor- 
rection was chosen on the basis of its prediction of a good 
null point in the 700 nm region. The detrital absorption 
spectra o^A) were then subtracted from the particulate 
absorption a p ( A) spectra to yield the phytoplankton absorp- 
tion spectra a ph ( A). 

[ 13 ] High-Performance Liquid Chromatography (HPLC) 
Pigments. Phytoplankton pigment concentrations were 
measured by HPLC using the procedure described by Van 
Heukelem and Thomas [2001], A total of 18 pigment 
groupings were identified in each of the samples (Table 2). 
Samples were processed in a manner similar to that for the 
absorption measurements and described earlier. 

2.2. Global Absorption Spectra and HPLC Pigment 
Data Set 

[ 14 ] An additional 608 phytoplankton absorption spectra 
and associated HPLC pigment measurements were obtained 
from the NASA SeaBASS data archive and added to the 
analysis data set [Werdell and Bailey, 2002; Werdell et ah, 
2003]. The SeaBASS archive contains over 13,692 absorp- 
tion spectra measurements and of these only 7,051 were 
determined usable for this study after assessing the spectral 
resolution, noise level, and other QA/QC concerns. Only 
608 of the total 7,051 available phytoplankton absorption 
measurements could be uniquely matched to one of the 
9,63 1 SeaBASS archive HPLC pigment measurements. An 
ensemble image, inclusive of the U.S. east coast data, of the 
various phytoplankton absorption spectra (Figure 3) shows 
that a wide range of spectral shapes is observed. Because of 
the large number of sources of these data it is not practical 
to provide a level of detail on the various laboratory pro- 
cessing involved in obtaining the phytoplankton absorption 
spectra and HPLC pigment concentrations. 

2.3. Determination of Individual Pigment Absorption 
Spectra and Pigments 

[ 15 ] Phytoplankton pigment composition is now routinely 
determined using standard HPLC analysis, which yields 
the concentrations of the majority of pigments. By com- 
bining these pigment concentrations with pigment-specific 
absorption spectra, it is possible to reconstruct the phyto- 
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Figure 3. Composite plot of the phytoplankton specific 
absorption spectra for all of the samples used in the analysis. 
The inset is a histogram of the log 10 of the specific absorp- 
tion at 440 nm. 

plankton absorption spectra for the algal culture/sample, 
such that 

m 

a 'ph(>) = y>,*(A), (4) 

i=i 

where c,- is the concentration of the individual pigments 
derived from HPLC analysis and af( A) is the pigment- 
specific (a.k.a. weight-specific) absorption coefficient for 
the phytoplankton pigment (Table 3). Bidigare et al. [1990] 
used this type of analysis in an attempt to reconstruct phy- 
toplankton absorption spectra using 5 dominant pigment- 
specific absorption spectra. This reconstruction technique 
can work well in open ocean regions where the package 
effect is low but typically does not work for areas such as 
coastal regions where the package effect is significant. 


[i6] When a large number of (n) related samples of phyto- 
plankton absorption spectra and HPLC observations are 
available it becomes possible to relate the pigment-specific 
absorption coefficients and HPLC pigment concentrations 
to the level of phytoplankton absorption measured at a spe- 
cific wavelength as 


/ c i=lj=l • • • Ci= m j= 1 \ 


f~aU (A) \ 


( a ph,j=l (A) ^ 

V Q= 1 J—n ' ' ' Ci =m j =n j 


\d* =m (A) / 


V &phj=n (^) / 


where c,- 7 - is the observed pigment concentration of the ith 
pigment and the jth sample, «f(A) is the derived pigment- 
specific absorption for the ith pigment, and a ph j(X) is the 
measured phytoplankton absorption for the jth sample and 
at a given wavelength (A). At this point the various con- 
centrations and absorption terms are members of a system of 
linear equations and can be rewritten in matrix form as, 

CA = A ph . (6) 

[ 17 ] We refer the solving the inverse of this equation to 
obtain an estimate of the unknown pigment-specific 
absorption coefficients for each of the phytoplankton pig- 
ments, a*(A), as “Ensemble Sample Deconvolution” [ESD] 
since it requires an ensemble of related phytoplankton 
absorption spectra and HPLC pigment measurements for 
proper implementation. The resulting pigment-specific 
absorption functions or spectra can then be used in combi- 
nation with the HPLC pigment concentrations to reconstruct 
the phytoplankton absorption spectra through the linear 
combination of the products of the individual pigment 
concentrations and pigment-specific absorption spectra as 
shown in equation (1). Care must be taken in order to avoid 
confusion between the pigment-specific absorption coeffi- 
cients obtained through laboratory measurements with those 
based on this inversion technique. They are not calculated in 
the same manner, one is based on laboratory measurements 
a*(A) and the other is a purely mathematical solution set, 
af(A). 

[is] In this paper, we compare the results for solving 
equation (6) using four numerical techniques. They include: 
Singular Value Decomposition (SVD) [Press et al ., 2007], 
Non-Negative Least squares (NNLS) [Lawson and Hanson, 


Table 3. Definition of Various Tenns 


Term 

Definition 

Units 

o p (A) 

Absorption coefficient of particles 

nT 1 

aA A) 

Absorption coefficient of methanol-extracted particles 

m 1 

A) 

Absorption coefficient of phytoplankton [a p ( A) — A)] 

nf 1 

^ phi A) 

Reconstructed absorption coefficient of phytoplankton 

nT 1 

a*,{ A) 

Pigment-specific absorption coefficient of phytoplankton pigment group i obtained from absorption 
measurement of in vitro pigment samples and corrected for shift in peaks 

m 2 mg -1 

a*,(A) 

Numerically estimated pigment-specific ‘absorption’ coefficient of phytoplankton group i obtained from 
Ensemble Sample Deconvolution (ESD) solutions. Not a true pigment- specific absorption term. 

m 2 mg -1 

a* p h( A) 

Absorption coefficient of phytoplankton normalized to the ‘unpackaged’ absorption estimates at 675 nm 

nT 1 

C; 

Concentration of phytoplankton pigment group i 

mg itT 3 

a so l A) 

Absorption coefficient of phytoplankton in the absence of the pigment package effect 

m _1 

®miss( A) 

Absorption coefficient of missing 


a 

Value of missing absorption coefficient at 675 nm 

m -1 

V 

Power function term in ‘missing’ absorption coefficient model 

n.d. 
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1974] and two Nonlinear Least Squared Minimization 
schemes (MRQMIN [Levenberg, 1944; Marquardt, 1963; 
Press et ciL, 2007] and LMDIF [Hiebert, 1980]). The utility 
of these numerical methods is that they are not impacted by 
potential spectral shifts of individual pigment absorption 
curves that does occur with direct absorption measure- 
ments of individual pigments dissolved in a solvent such as 
methanol, as pointed out by Bidigare et al. [1990] and 
Bricaud et al. [2004]. In addition, these methods estimate 
the in vivo rather than in vitro pigment-specific absorption 
spectra, thereby explicitly including pigment package effects. 

[ 19 ] Singular Valued Decomposition (SVD). The first of 
these methods is approached by taking the inverse of C, 
thereby making it possible to directly solve the matrix 
equation in order to determine unique absorption spectra for 
individual pigments such that 

C 1 (CA) = C-'Aj,,, = I„A = A. (7) 

In practice, because of the difficulties of direct matrix 
inversion and problems related to singular matrices, the 
approach to this method is often to use the Singular Value 
Decomposition (SVD) methodology in order to solve for the 
unknown pigment-specific absorption values. A thorough 
presentation of this method is found in the work by Press 
et al. [2007]. 

[ 20 ] The pitfalls that can occur in the SVD approach is 
that the solution set cannot be guaranteed to be completely 
positive, so that negative ‘absorption' values can result. One 
needs to approach the results from such analysis with an 
understanding that these negative values are not related to 
something akin to fluorescence but are related to the solu- 
tion noting a negative influence from that specific pigment 
in relation to the measured phytoplankton absorption at that 
wavelength. While the units of the solutions are in nr mg -1 , 
the values must be thought of as a component to the phy- 
toplankton absorption value and not a separate process. 

[ 21 ] Non-Negative Least squares (NNLS). In order to 
address the issue of possible negative absorption coefficients 
within the solution, a second method called Non-Negative 
Least squares [Lawson and Hanson , 1974] is carried out. In 
this approach, the previous linear system of equations is 
reformulated to 

minimize || CA — A p i, || subject to A> 0. (8) 

The details to this technique are presented by Lawson and 
Hanson [1974], The essence of the technique is that a 
solution is found for the various unknown pigment-specific 
absorption coefficients, a* (A) in a manner similar to the 
SVD technique but with the caveat that all the coefficients 
are positive. 

[ 22 ] Nonlinear Least-Squared Minimization (MRQMIN/ 
LMDIF). A third strategy that is used to obtain estimates of 
the pigment-specific absorption spectra is to apply a non- 
linear minimization algorithm to the Sum of Squared Error 
(SSE) between the total absorption and the estimated total 
absorption such that the 

n / m 

SSE = I a p i,j( A) — ^2 c ij 
j= 1 \ '=1 


By including a large number of field samples (n) that 
have been analyzed for absorption spectra and the suite 
of HPLC pigments, it is possible to develop the pigment- 
specific absorption spectra through the solution of these 
matrix equations at each wavelength. Taking the absolute 
values of the pigment-specific absorption coefficients forces 
the solutions to be positive. Two different algorithms are 
implemented. This nonlinear minimization scheme is called 
Levenberg-Marquardt [ Levenberg , 1944; Marquardt, 1963] 
and is implemented using both the MRQMIN algorithm 
described by Press et al. [2007] and the LMDIF algorithm 
developed at Argonne National Laboratory [Hiebert, 1980], 

[ 23 ] The complete set of 874 phytoplankton absorption 
spectra and pigment HPLC measurements are used to derive 
specific absorption coefficients using the three (SVD, 
NNLS, MRQM1N/LMD1F) methods outlined above. One 
set was created using solely the 266 samples from the U.S. 
east coast data set and the other from a combination of the 
full 874 observations, allowing for a comparison of regional 
versus global solutions. These are presented in the results 
section. The purpose of using multiple inverse model ap- 
plications (SVD, NNLS, MRQM1N/LMDIF) is to verify 
that the inverse model solutions are robust and not influ- 
enced by the choice of matrix inversion application. Vallino 
[2000] used 12 different minimization schemes to optimize 
a system of coupled ordinary differential equations that 
simulated an ocean ecosystem. All of the minimization 
schemes converged on different model parameter sets. 

2.4. Estimation of HPLC Pigment Concentrations 

[ 24 ] Once estimates for pigment-specific absorption 
coefficients are available, either through laboratory mea- 
surements [Bidigare et al., 1990; Bricaud et al., 2004] or 
through numerical analysis such as the ESD approaches 
outlined above, it becomes possible to combine those with 
phytoplankton absorption spectra in order to estimate the 
HPLC pigment concentrations. By expanding upon the 
phytoplankton absorption spectra reconstruction technique 
of Bidigare et al. [1990], the phytoplankton absorption can 
be written as 


/at i(A=l) ... af =M ( A = 1) \ 


( ° i=X \ 


( — 1) ^ 

\a? = i(A = L) ••• a* = hi{A = L) ) 


V Cj=M J 


V &ph(^ / 


( 10 ) 

where df(A) is the estimated pigment-specific absorption of 
the ith pigment at a specific measured wavelength, Cj is the 
estimated concentration of pigment j, and a p /,( A) is the total 
measured absorption due to phytoplankton. At this point the 
equations are a system of linear equations and can be 
rewritten in matrix form as, 

AC = A ph , (11) 

where the matrix A is the set of pigment-specific absorption 
coefficients, C is the array of pigment concentrations to be 
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Pearson Product Moment Coefficient of Correlation [r] 
(a) 0.000 I 1.000 



Pigment Name 


Pearson Product Moment Coefficient of Correlation [r] 
(b) O.OOOH [ 1.000 



Figure 4. Plot showing the Pearson Product Moment Coefficient of Correlation. The names of the var- 
ious pigment types are shown along the bottom and left of the plot. The color bar scale is linear. Truncated 
correlation values are written within the associated color block and those correlation values that are not 
significantly greater than zero are left blank. 


estimated and A ph is the phytoplankton absorption spectra. At 
this point the equation looks similar to that of equation (6) 
and so similar inversion techniques can be applied to esti- 
mate the HPLC pigment concentrations using phytoplankton 
absorption spectra measurements. 

[ 25 ] In practice, independent data sets were created by 
randomly splitting the U.S. east coast (266 samples) and 
global (874 samples), into two halves each. Phytoplankton 
absorption spectra and HPLC observations from one of the 
splits was used to calculate the pigment-specific absorption 
spectra. The results from this inversion were used with the 
phytoplankton absorption spectra from the second split data 
set to obtain an independent estimate of the phytoplankton 
pigment suites which were then compared to the HPLC 
observations from the second split. 

3. Results 

3.1. HPLC Pigment Correlation 

[ 26 ] The Pearson Product Moment Coefficient of Corre- 
lation, r, is calculated for each set (U.S. east coast and 
combined global) of HPLC pigment observations (Figure 4). 
Correlations between pigments from the U.S. east coast data 
set (Figure 4a) were much higher overall than that observed 
within the larger global HPLC data set (Figure 4b). High 
levels of correlation between pigments is a concern because 
it can cause multiple inverse model solutions. 

3.2. Pigment-Specific Absorption Spectra/Coefficients 

[ 27 ] Four different solution sets were obtained from car- 
rying out the inversion to obtain pigment-specific absorp- 
tion spectra, one for each of the numerical methods applied 
(SVD, NNLS, MRQMIN, LMD1F). Solutions for the SVD 
and NNLS are shown in Figures 5a and 5b, respectively. In 


the SVD application, negative coefficients are possible 
because the method is a pure matrix inversion. The other 
three applications (NNLS, MRQMIN, and LMDIF) are 
nonlinear least square applications that allows the solutions 
to be constrained to be positive. What is obvious when 
comparing between the four solutions is that: i) the SVD 
solution contains negative values; ii) only the SVD and 
NNLS solutions vary smoothly across the spectra; iii) the 
NNLS, MRQMIN and LMDIF solutions show similar 
spectral shapes with carotenoid having the largest specific- 
absorption levels overall; and, iv) the carotenoid specific- 
absorption is almost an order of magnitude higher than any 
of the laboratory measured in vitro specific-absorption 
spectra (Figure 1). In addition to the variability between the 
solution sets, it needs to be pointed out that the solutions are 
obtained by fitting ‘packaged’ absorption spectra to the 
‘unpacked’ HPLC pigment measurements. Because HPLC 
pigment observations show high levels of correlations 
between various pigments and because of the wide range in 
possible levels of pigment package effect there should be no 
expectations that the various pigment absorption spectra 
model solutions would reflect those obtained from labora- 
tory measurements, such as shown in Figure 1. 

[ 28 ] The various solutions to the pigment-specific 
absorption spectra were used with the observed HPLC 
pigment measurements in equation (1) to reconstruct the 
phytoplankton absorption spectra in order to provide for a 
way to quantitatively assess the ability of these inverse 
model solutions (Figure 6). In addition to these solutions, 
a ‘baseline’ solution was established using the pigment- 
specific absorption spectra from Bricaud et al. [2004]. The 
coefficient of determination (r 2 ) for each of these solutions 
as a function of wavelength is shown in Figure 7. The 
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solutions, though seemingly limited due to its rendering of 
negative absorption coefficients, outperform the three other 
inversion algorithms (NNLS, MRQMIN, and LMDIF). The 
three least squares inversion algorithms show very similar 
linear regressions slopes and confidence intervals (Figure 7a) 
and coefficients of determination (Figure 7b). All solutions, 
including the SVD and ‘baseline’ reconstruction show a 
marked decrease in r 2 values in the range of 550 to 600 nm. The 
correlations and nonnalized standard deviations of the various 
inverse model solutions (Figure 8) shows that the correlations 
for all of the solutions is high (>0.9) and that the standard 
deviations are slightly lower than the observations. 

[ 29 ] Much of the present difficulty of successful spectral 
inversions is due to the influence of the pigment package 
effect on phytoplankton absorption spectra. This phenome- 
non can be observed by comparing the absolute absorption 
measured at 675 nm with the concentration of chlorophyll a, 
which is the primary absorbing pigment at that wavelength 
(Figure 9). The falloff from the ‘expected’ reconstructed 
absorption levels (shown as the black line in Figure 9 and 
having a slope equal to the chlorophyll a-specific absorption 



NNLS Solution 

_i I l I I I I I I I — I I I I I — I I I I I I I I I 
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Fucoxanthin 

Neoxanthin 

Violaxanthin 

19 Hexa...fuco. 
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Figure 5. Pigment-specific absorption spectra/modes for 
the 18 pigments analyzed in this study obtained from the 
(a) SVD, and (b) NNLS algorithms. Solutions for the 
MRQMIN and LMDIF algorithms (not shown) are similar 
to that of the NNLS solution. 


baseline solution using the Bricaud et al. [2004] pigment- 
specific absorption spectra is limited to the spectral range 
of 400 to 700 nm, whereas reconstructions from the four 
inversions can extend into the UV region. All inverse model 
solutions show a dramatic increase in the ability to recon- 
struct the total phytoplankton absorption spectra relative to 
the baseline Bricaud et al. [2004] reconstructions. The SVD 


Figure 6. A comparison of the absorption spectra observa- 
tions (red curve) and six absorption spectra reconstructions 
for one of the 792 absorption spectra and pigment observation 
pairs. The SVD, MRQMIN, LMDIF and NNLS solutions 
(thin green, teal, light and dark blue curves, respectively) 
extend across the full 300 to 700 nm spectral range. The dif- 
ference (thick black curve) between the observations (thin red 
curve) and the Bricaud et al. [2004] reconstruction (thin black 
curve ranging from 400 to 700 nm) appears to look like a 
power function. This difference is the ‘missing absorption 
term’ that was modeled using a power function (thick mauve 
curve). The sum (thick purple curve) of the standard Bricaud 
et al. [2004] (thin black curve) reconstruction and the fitted 
power function (thick black curve) compares very well to the 
observations. 
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at 675 nm) is primarily a result of the pigment package effect. 
In order to assess this, all the measured phytoplankton 
absorption were nonnalized to 675 nm by multiplying the 
absorption spectra by a normalization term that can be 
thought of as the inverse pigment package effect term at 
675 nm, such that the spectra unpackaged, relative to 675 nm, 
or ‘normalized’ absorption spectra, 


^J2 c ia *( 675 nm)'' 


Gphity Uph (A) 


V 


Up/, (675 nm) 


u^,/,(A) 


a '( 675 nm) 
a p h[615 nm) 

(12) 



Wavelength [nm] 


(b) 



Wavelength [nm] 


[ 30 ] Estimates for a'{ A = 675 nm) were obtained using 
pigment-specific absorption values from Bricaud et al. 
[2004], Once all of the absorption spectra were ‘normal- 
ized’ (Figure 10) the various inversion algorithms described 
above (SVD, NNLS, MRQMIN, LMDIF) were reapplied to 
these ‘normalized’ absorption spectra. The statistical anal- 
ysis on the resulting absorption spectra (Figures 7a and 7b) 
shows a dramatic improvement as compared to those obtained 
from inversions done using the original measured absorption 
spectra. Using this as a guide, the Bidigare etal. [1990]model 
was reformulated to include this ‘normalization’ term 
(equation (12)) such that the modified the reconstructed 
absorption spectra, 


a P h( x ) 


J* P h(615 nm) 
f)c,a,*( 675 nm) '=1 


(13) 


The spectral reconstructions carried out by Bricaud et al. 
[2004] systematically showed lower absorption than the 
measured absorption spectra, leading them to argue that a 
missing absorption material was present in the phytoplankton 
samples. This is also true for many of the absorption spectra in 
our data set. In the majority of the samples from this study, 
the spectral reconstructions following Bricaud et al. [2004] 
showed higher absorption than the measured absorption 
spectra. This difference is expected since the Bricaud et al. 
[2004] samples were primarily obtained in open ocean set- 
tings, where the pigment package effect is minimal and a 
majority of the samples analyzed in this present study are 
from coastal ocean areas where the pigment package effect is 
often very significant. However, upon reconstruction of the 
absorption spectra using equation (13), which is equivalent to 
‘normalization’ of the absorption spectra to 675 nm (Figure 
10), the spectral reconstructions show lower absorption than 
the measured absorption spectra. 


Figure 7. The (a) slopes and confidence intervals, and 
(b) coefficient of determination (r 2 ) from comparison of the 
reconstructed versus measured phytoplankton absorption spec- 
tra for all samples. The ‘baseline’ level from using the Bricaud 
et al. [2004] absorption spectra (black curve) is limited to 
between 400 to 700 nm. The SVD solutions (red curves) are 
shown for solutions to the reconstructions using the raw 
(solid red curve) versus spectra nonnalized to 675 nm (dashed 
red curve). The NNLS, MRQMIN and LMDIF solutions 
(blue curves) are shown for solutions to the reconstructions 
using the raw (solid blue curves) versus the spectra nonnalized 
to 675 nm (dashed blue curves). The r 2 values of spectral 
reconstructions using the Bricaud et al. [2004] to total phyto- 
plankton absorption spectra ‘normalized’ to 675nm (dashed 
black curves) also show a dramatic improvement over the 
‘baseline’ reconstructions (solid black curves). Further 
improvements to this modified reconstruction are possible by 
fitting the coefficients to a logarithmic model for the observed 
missing absorption component (solid fuchsia curves). 
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Figure 8. Taylor plot of the reconstructed phytoplankton 
absorption spectra normalized to the observations. The color 
of the dots indicates the type of inverse model technique 
used. See Tables 4 or 5 for a listing. 


[ 31 ] Inspection of the difference between the ‘nonnalized’ 
measured absorption and the reconstruction spectra shows 
that the gross spectral shape of this missing absorption 
resembles a logarithmic function, 

( 675 nm\ ^ 

dm«i(A) = a( — ^ — J . (14) 


In order to diagnose this further, the original Bidigare et al. 
[1990] total absorption spectra reconstruction equation 
(equation (1)) was further modified to include these new 
terms such that 


a'phW = ' 


a ph ( 675 nm) 


a + c i a f (675 nm) 

i=i 


( 675 nm\ v * 

+2>a,*(A) 


(15) 


[ 32 ] This equation was used with the observations of 
phytoplankton absorption and pigment concentration to fit 
the two parameters (a, ip) in the missing absorption term, 
resulting in a dramatic increase (r 2 > 0.98) in the ability to 
reconstruct phytoplankton absorption spectra (Figures 7a 
and 7b). In addition, one coefficient, a, shows strong cor- 
relation with the chlorophyll a (Figure 1 1 ) while the slope of 
the logarithm, ip, (not shown) did not. 

[ 33 ] Estimation of the pigment package effect at 675 nm 
has previously been done by Nelson et al. [1993] by using 
the quotient of the reconstructed versus the measured 


absorption values at 675 nm and this method of estimation 
has been further argued by Bricaud et al. [2004] which notes 
that the “package effect index,” 

&*( A)=^§, (16) 

where a ptl ( A) is the actual absorption coefficient and a so/ ( A) 
is the absorption coefficient from the same material when 
dispersed within a solution [Morel and Bricaud, 1981]. By 
making use of the fully reconstructed solutions through the 
modification of equation (16) to remove the ‘normalization’ 
term, the pigment package effect can be estimated as, 

n*f'\ _ a phW _ a ph (A) / 17 \ 

Ua I'D — — m l 1 'J 

asol[\) Q (675_nm^ + ^ c . a *(A) 

;=i 

for the entire visible spectrum (Figure 12). When the 
missing absorption components are accounted for in the 
absorption spectra, the estimated pigment package effect 
term shows little spectral variability and ranges between 
samples across its full range of 0 to 1, with a mean near 0.5, 
such that for any one spectra 

Qa( A)~ £ ~ const. (18) 

a + Y2 Cjdf (675 nm) 

1=1 

When the pigment package effect is calculated using 
equation (16) and the initial reconstructed phytoplankton 
absorption spectra the values are often higher than 1 and 
show much more variability across the spectra. 



Figure 9. Observed phytoplankton absorption levels at 
675 nm versus chlorophyll a concentrations. The line denotes 
the expected level of absorption from unpackaged chloro- 
phyll a absorption. 
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Figure 10. Composite plot of the phytoplankton absolute 
absorption spectra normalized to 675 nm for all of the sam- 
ples used in the analysis. 


3.3. Inversion of the Pigment-Specific Absorption 
Spectra to Estimate Phytoplankton Pigments 

[ 34 ] A goal in modeling phytoplankton absorption spectra 
is to make use of the reconstruction models to estimate 
phytoplankton pigments directly from the phytoplankton 
absorption spectra (equation (11)). The same numerical 


inversion techniques used to solve for the pigment-specific 
absorption spectra can also be used to solve these sets of 
equations. A set of five inverse model solutions was carried 
out on both the raw and ‘normalized’ total absorption spectra. 
These include: a) SVD, b) SVD-NNLS (which solved the 
equations using the NNLS technique on the smooth pigment- 
specific absorption spectra solutions from the SVD inversion, 
c) NNLS, d) MQRMIN, and e) LMD1F. The coefficients of 
determination (r 2 ) are shown in Table 4 for all solutions using 
the U.S. east coast (COBY+BIOME+MAA) observations 
and in Table 5 for all solutions using the combination of all 
observations (SeaBASS+COBY+BIOME+MAA). The pig- 
ment-specific absorption spectra required for the inverse 
models was obtained from the inverse model solutions for 
these spectra using a randomly split portion of the full data set. 
The measured absorption spectra and HPLC pigments from 
the other part of the data was used to independently validate 
the inverse model solutions. 

[ 35 ] U.S. East Coast Pigment Estimates: In all cases, the 
concentrations being estimated reflect those measured in the 
HPLC pigment analysis. Chlorophyll a, carotenoids, fiico- 
xanthin, violaxanthin, diadinoxanthin and peridinin all show 
some level of adequate correlations (r > 0.6), though 
fucoxanthin and the carotenoids are the only two pigments 
that are retrieved consistently across all the inversion 
methods. Chlorophyllide, alloxanthin, diatoxanthin, zeax- 
anthin, lutein and chlorophyll b all have poor r 2 values. 
Finally, both 19'-butanoyloxyfucoxanthin (but-fuco) and 
19'-hexanoyloxyfucoxanthin (hex-fuco), two poorly corre- 
lated pigments in the HPLC measurements (Figure 4) also 
had very poor retrievals across all inversions. Inversion 
solutions that used the 675 nm normalized spectra had much 
better predictions of pigments, with chlorophyll a values 
being almost exact. 



Figure 11. Magnitude of the missing absorption term, a, 
versus the observed HPLC chlorophyll a concentrations. 
This term is the absolute level of missing absorption at 
675 nm. 



400. 450. 500. 550. 600. 650. 700. 


Wavelength [nm] 

Figure 12. Estimates of the pigment package effect for all 
samples calculated using the ratio of measured in vivo 
absorption to modeled in vitro absorption spectra. 
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Table 4. Coefficient of Determination (r 2 ) Levels for Retrieval of Pigment Concentrations Using Only the COBY, BIOME and MAA 
Data Sets'* 


Pigment 

Name 



Raw 



Auto-Normalized b 


Normalized to 675 nm 


SVD 

SVD-NNLS 

NNLS 

MRQMIN 

LMDIF 

Bricaud 

Bricaud + 

Missing SVD SVD-NNLS 

NNLS 

MRQMIN 

LMDIF 

chlorophyll c (cl+c2+c2) 

0.56 

0.69 

0.27 

0.38 

0.40 

0.89 

0.57 

0.78 

0.86 

0.48 

0.38 

0.47 

chlorophyllide 

0.36 

0.47 

0.23 

0.47 

0.05 

(b) 

(b) 

0.42 

0.45 

0.05 

0.00 

0.01 

pheophoribide 

0.51 

0.59 

0.22 

0.24 

0.23 

(b) 

(b) 

0.45 

0.58 

0.31 

0.49 

0.25 

Peridinin 

0.60 

0.65 

0.54 

0.52 

0.47 

0.16 

0.25 

0.48 

0.60 

0.30 

0.29 

0.36 

19’ butanoy-loxyfucoxanthin 0.31 

0.42 

0.53 

0.46 

0.56 

0.01 

0.10 

0.39 

0.52 

0.50 

0.30 

0.69 

fucoxanthin 

0.84 

0.87 

0.78 

0.80 

0.80 

0.94 

0.70 

0.89 

0.93 

0.58 

0.65 

0.59 

neoxanthin 

0.36 

0.41 

0.27 

0.43 

0.37 

(b) 

(b) 

0.45 

0.57 

0.33 

0.36 

0.39 

violaxanthin 

0.63 

0.72 

0.53 

0.63 

0.57 

(b) 

(b) 

0.47 

0.74 

0.43 

0.07 

0.29 

19' hexanoy-loxy fucoxanthin 0.30 

0.40 

0.00 

0.40 

0.17 

0.31 

0.14 

0.37 

0.52 

0.01 

0.13 

0.15 

diadinoxanthin 

0.71 

0.71 

0.37 

0.41 

0.42 

0.74 

0.51 

0.73 

0.74 

0.42 

0.29 

0.24 

alloxanthin 

0.31 

0.41 

0.11 

0.01 

0.20 

0.50 

0.30 

0.48 

0.58 

0.13 

0.18 

0.12 

diatoxanthin 

0.15 

0.26 

0.00 

0.01 

0.01 

(b) 

(b) 

0.25 

0.45 

0.32 

0.05 

0.30 

zeaxanthin 

0.36 

0.44 

0.18 

0.22 

0.01 

0.06 

0.20 

0.27 

0.34 

0.04 

0.00 

0.01 

lutein 

0.20 

0.26 

0.27 

0.18 

0.08 

(b) 

(b) 

0.18 

0.18 

0.07 

0.24 

0.32 

chlorophyll b c 

0.27 

0.33 

0.23 

0.46 

0.38 

0.00 

0.15 

0.28 

0.53 

0.41 

0.34 

0.42 

chlorophyll a c 

0.76 

0.81 

0.59 

0.36 

0.39 

NA d 

NA d 

1.00 

1.00 

0.98 

0.86 

0.96 

phaeophytin a 

0.59 

0.65 

0.53 

0.45 

0.48 

(b) 

(b) 

0.63 

0.72 

0.37 

0.48 

0.42 

carotenoids 

0.74 

0.83 

0.78 

0.83 

0.79 

0.57 

0.00 

0.92 

0.94 

0.68 

0.83 

0.78 


“Values greater than 0.5 are in bold. 

b Entries with (b) indicate pigment not included in inversion due to lack of pigment-specific absorption coefficients. 
“Includes both mono and divinyl. 

d Inversion solutions normalized against chlorophyll a values, making correlations irrelevant. 


[36] Because of the very good agreements from fitting the 
observed absorption spectra with the modified Bidigare et al. 
[1990] absorption spectra model (equation (15) and Figure 7), 
the capability of inverting this new model estimate to the 
HPLC pigment observations was also assessed. For com- 
parison, an additional inversion was carried out using 
equation (13), which does not contain the missing absorption 
term. Unfortunately, because of the way these equations are 
constructed it is not possible to obtain quantitative estimates 
due to the self-scaling nature of the pigment package effect 


term, any changes in the estimated pigment concentrations 
are readily offset by a balanced change in the packing term, 
yielding similar results. As a consequence, without additional 
information on either the level of pigment package effect or 
the actual concentration of a key pigment, only relative pig- 
ment ratios can be retrieved. Therefore, in order to yield some 
quantitative assessment of using this equation the chlorophyll 
a concentration was set to the HPLC observations and the 
LMDIF algorithm was used to solve for the remaining 
unknown pigment concentrations. In both cases (Table 4, 


Table 5. Coefficient of Detennination (r°) Levels for Retrieval of Pigment Concentrations Using Randomly Split Data Set of All the 
Absorption Spectra and HPLC Observations'* 


Raw Auto-Normalized b Normalized to 675 rim 


Name SVD SVD-NNLS NNLS MRQMIN LMDIF Bricaud Bricaud + Missing SVD SVD-NNLS NNLS MRQMIN LMDIF 


chlorophyll c (cl+c2+c2) 

0.01 

0.00 

0.00 

0.00 

0.01 

0.17 


0.10 

0.23 

0.05 

0.15 

0.06 

chlorophyllide 

0.03 

0.12 

0.03 

0.00 

0.04 

(b) 

(b) 

0.02 

0.26 

0.06 

0.10 

0.09 

pheophoribide 

0.02 

0.02 

0.03 

0.02 

0.02 

(b) 

(b) 

0.01 

0.06 

0.09 

0.14 

0.15 

peridinin 

0.06 

0.16 

0.08 

0.16 

0.16 

0.54 


0.03 

0.32 

0.22 

0.26 

0.28 

19' butanoy-loxyfucoxanthin 

0.01 

0.00 

0.01 

0.01 

0.02 

0.00 


0.00 

0.01 

0.07 

0.00 

0.03 

fucoxanthin 

0.15 

0.18 

0.14 

0.29 

0.33 

0.89 


0.31 

0.70 

0.49 

0.66 

0.66 

neoxanthin 

0.01 

0.16 

0.05 

0.27 

0.10 

(b) 

(b) 

0.00 

0.06 

0.01 

0.03 

0.04 

violaxanthin 

0.16 

0.30 

0.00 

0.55 

0.03 

(b) 

(b) 

0.07 

0.30 

0.43 

0.50 

0.25 

19' hexanoy-loxy fucoxanthin 

0.00 

0.10 

0.07 

0.07 

0.04 

0.00 


0.04 

0.16 

0.04 

0.15 

0.06 

diadinoxanthin 

0.04 

0.07 

0.00 

0.15 

0.08 

0.35 


0.00 

0.43 

0.31 

0.58 

0.42 

alloxanthin 

0.24 

0.46 

0.53 

0.44 

0.39 

0.18 


0.25 

0.52 

0.34 

0.52 

0.57 

diatoxanthin 

0.30 

0.53 

0.01 

0.00 

0.41 

(b) 

(b) 

0.08 

0.33 

0.08 

0.03 

0.17 

zeaxanthin 

0.39 

0.63 

0.44 

0.72 

0.65 

0.37 


0.22 

0.56 

0.76 

0.77 

0.71 

lutein 

0.01 

0.11 

0.22 

0.26 

0.19 

(b) 

(b) 

0.00 

0.03 

0.07 

0.09 

0.09 

chlorophyll b c 

0.03 

0.39 

0.39 

0.60 

0.30 

0.24 


0.09 

0.50 

0.22 

0.54 

0.34 

chlorophyll a c 

0.22 

0.21 

0.44 

0.50 

0.52 

NA d 

NA d 

0.81 

1.00 

0.993 

0.998 

0.997 

phaeophytin a 

0.26 

0.31 

0.22 

0.20 

0.15 

(b) 

(b) 

0.04 

0.16 

0.13 

0.14 

0.15 

carotenoid 

0.02 

0.08 

0.02 

0.05 

0.06 

0.09 


0.23 

0.33 

0.39 

0.49 

0.46 


“Values greater than 0.5 are in bold. 

b Entries with (b) indicate pigment not included in inversion due to lack of pigment-specific absorption coefficients. 
“Includes both mono and divinyl. 

d Inversion solutions normalized against chlorophyll a values, making correlations irrelevant. 
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Auto-Normalized cases), the ability to estimate the measured 
HPLC pigments was worse than what was encountered in all 
of the model inversions that used results from the simple 
matrix inversions to obtain pigment-specific absorption 
spectra or more appropriately “absorption modes.” 

[ 37 ] Global Pigment Estimates: Compared to the U.S. east 
coast pigment estimates, using a more globally distributed 
data set yielded coefficients of determination (r 2 ) that were 
significantly lower (Table 5). Carotenoid estimates showed 
the most significant drop in performance. Chlorophylls a 
and b, zeazanthin, alloxanthin, violaxanthin, and focux- 
anthin had reasonable (>0.5) estimations but did not perform 
consistently among the various inverse models. Chlorophyll 
a predictions were again very strong (>.8) in the inverse 
model solutions that used the normalized absorption spectra 
as the inputs to the inverse models. 

4. Discussion 

[38] Developing models that link phytoplankton pigments 
to the lOPs of phytoplankton absorption spectra is an 
essential step for advancing the capabilities of hyperspectral 
inverse models for ocean remote sensing reflectance 
observations. The inversions above show that it is possible 
to extract pigment-specific absorption ‘modes,’ that can be 
used with measured pigment concentrations to reconstruct 
observed phytoplankton absorption spectra. Unlike the 
presently available laboratory derived pigment-specific 
absorption spectra [Bidigare et al., 1990; Bricaud et al, 
2004], which are limited in spectral range between 400 to 
700 nm, the ESD technique is able to derive absorption 
relationships/modes well into to ultraviolet (UV) region of 
the spectrum. However, the coefficients of determination 
(r 2 ) for this spectral region (300-400 nm) is much less 
(ranging from about 0.5 to 0.92, not shown) than for the 
visible region of the spectrum (Figure 7). The reason for 
this decrease is likely due to the fact that the pigments that 
were used in the inversion process did not contain the 
key pigments that absorb strongly in the UV region, such 
as mycosporine-like amino acids (MAAs) [Moisan et al. , 
2010 ]. 

[ 39 ] Reconstructing Phytoplankton Absorption Spectra. 
The differences between the ESD reconstructions and the 
Bricaud et al. [2004] reconstruction that was modified to 
include a power function to resolve the likely missing 
absorption term (equation (15)) presents some puzzling 
observations. The first of these is related to whether the 
pigment package effect was adequately estimated by the 
inclusion of a ‘normalization’ term in equation (15). The pre- 
vious work of Bricaud et al. [2004] argues for a missing 
absorption component. This missing absorption compo- 
nent was also observed in phytoplankton absorption spectra 
collected at the depth of the 1 % light level off the coast of 
California [ Nelson et al, 1993], Certainly the lack of: a) a 
complete set of pigment-specific absorption spectra (only 
14 versus the 18 measured HPLC pigments); b) MAA mea- 
surements; and, c) pigment concentrations and absorption 
spectra for the water soluble phycobiliprotiens (phycocyanin 
and phycoerythrin), which are absent from HPLC mea- 
surements, contributed to lowering the capabilities of all 
spectral reconstruction techniques. MAAs primarily absorb 
in the UV region, phycocyanin has peak absorption in the 


608-620 nm region [Glazer et al., 1973; Debreczeny et al., 
1993] and phycoerythrin absorbs primarily in the 550 to 
650 nm region [Bryant, 1982; Ong et al., 1984], Those 
spectral regions yielded lower r 2 values for all absorption 
spectra reconstructions. The missing absorption component 
might result from the presence of intracellular absorbing 
materials that are not HPLC measured pigments and the 
incomplete removal of the ‘detrital’ absorption components, 
which generally have featureless absorption spectra increas- 
ing to shorter wavelengths. Results from Bidigare et al. 
[1989] noted that these ‘non-photosynthetic chromophores’ 
affected quantum yield measurements that relied on phyto- 
plankton absorption spectra reconstruction techniques. The 
measured spectra might also be influenced by a spectrally 
varying pigment package effect that would alter the overall 
shape of the entire spectrum, thereby altering the shape and 
magnitude of the ‘missing 1 absorption component. In addi- 
tion to those issues mentioned above, there are also errors 
associated with the actual absorption spectra measurements 
[Mitchell et al., 2000], HPLC measurements [Bidigare et al, 
2005] and pigment extraction efficiencies. 

[ 40 ] A comparison (not shown) of all of the power func- 
tions against the difference between the Bricaud et al 
[2004] reconstructions and the observed absorption spectra 
(e.g., Figure 6) shows that while the power function does 
a fairly decent job at modeling the overall shape of the 
differences, the spectral absorption in the region between 
450 to 500 nm was consistently underestimated and over- 
estimated in the region between 600 to 625 nm. This con- 
sistent bias shows up in the linear regressions between the 
measured and reconstructed phytoplankton absorption 
spectra (Figure 7a). This suggests again that other pigments 
and their absorption spectra need to be included in the 
reconstructions or that the observed absorption spectra 
shapes are being modified by the spectral variability of the 
pigment package effect. The difference spectra in the Bricaud 
et al. [2004] reconstructions demonstrated carotenoid-like 
spectral signatures, which support our overall analysis. A 
power function fitting such spectrum would underestimate 
the absorption in the region from 450 to 500 nm. 

[ 41 ] It is interesting to note also that the inversion 
reconstructions within the UV region (300 to 400 nm, not 
shown) are actually fairly well represented, even though no 
information is available about MAA concentrations or their 
absorption spectra. Some of this is likely due to the fact that 
MAAs have been shown to be correlated to levels of the 
pigment diatoxanthin [Moisan et al, 2010], which is 
represented in the HPLC data set. It is likely that these 
various levels of co-variability between the pigment con- 
centrations, either observed (Figure 4a) or not (as in the 
MAAs), play a significant role in determining the final 
pigment-specific absorption spectra. 

[ 42 ] A noticeable decline in performance of all absorption 
spectra reconstructions occurred in the spectral region 
between 550 and 650 nm. This is likely caused by not 
including the contributions of phycobilin pigments to the 
phytoplankton absorption spectra reconstructions. Phycobi- 
lins were not included in this study because they are not part 
of the standard pigment suite measured in HPLC analysis. 
Finally, there is a consistent drop in r 2 near 700 nm for all 
spectral reconstruction methods attempted in this study. 
Certainly in this region the absorption levels are very low 
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compared to most other parts of the absorption spectra, 
which can enhance errors. 

[ 43 ] Pigment Package Effect. The pigment package effect 
spectra have traditionally been calculated by dividing 
observed phytoplankton absorption spectra by the reconstructed 
absorption spectra [. Nelson et al., 1993; Bricaud et al., 2004]. 
This is not possible for samples that demonstrate higher 
values of absorption that the reconstructed absorption spectra 
would suggest, yielding pigment package effects with values 
higher than 1, values outside the defined 0-1 range. Our 
results suggest that the missing absorption components 
should be part of the reconstructed absorption spectra equa- 
tion. In doing so, the spectral shape of the pigment package 
effect (Figure 12) for all of the samples shows modest vari- 
ability across the spectrum but a wide variability between 
the samples. In the grossest sense, reduced (higher value) 
levels of pigment package effect occurred around 480 and 
680 nrn, with higher levels near the 580 nrn region of the 
spectrum. This variability points back to the issue of not 
completely reconstructing the absorption spectra due to the 
variety of reasons listed previously. The manner in which the 
pigment package spectra are calculated (equation (18)) also 
causes the spectra to contain the measure of all of the errors 
from the absorption spectra reconstructions, making it much 
more difficult to believe the result. Geiderand Osborne [1987] 
used an alternate method for estimating the magnitude of 
the pigment package effect that entailed comparing the ratio 
of intact to disrupted cell absorption. Their results showed 
a range in values from 0.5 at the 435 nm, 1 .0 at 600 nm, and 0.7 
at 670 nm, which is markedly different from the general 
spectral variations observed in this study. Other examples of 
the pigment package effect estimated by the ratio of intact 
to disrupted cell absorption show a wide range of spectral 
variability [Osborne and Geider, 1989; Kirk, 1994]. One 
interesting example of a comparison between the absorption 
spectra of intact versus disrupted cells of a Synechococcus 
sp. (WH7803 clone) shows that the spectral shape of the 
pigment package effect be similar is shape to that of the 
‘missing absorption 1 term placed into the absorption model in 
equation (18) [see Osborne and Geider, 1989, Figure 2A], 
In the methodology that is used in this paper, an assumption 
is made in order to obtain estimates of the pigment package 
effect. The assumption is that difference between the 
‘unpackaged’ spectra and the standard Bricaud et al. [2004] 
reconstruction was characteristic of the absorption due to 
a missing absorption component that could be modeled fairly 
well using a power function that, oddly enough, is similar to 
a scattering function. However, given that MAA, phycocyanin 
and phycoerythrin absorption components were not included 
in the spectral reconstruction, and there have been a number 
of arguments for missing absorption [Nelson et al., 1993; 
Bricaud et al. 2004] and scattering [Osborne and Geider, 
1989] components, much work on this issue remains. Is this 
‘missing absorption term’ noted previously by Bricaud et al. 
[2004] and observed again in the results above due to 
absorption, some component of the pigment package effect, 
scattering, or a combination of them all? And, the issue of 
errors within the pigment absorption reconstructions have a 
direct impact on the spectral variability of the pigment package 
effect. 


[ 44 ] Estimating Phytoplankton Pigments using Absorption 
Spectra. The ability to reconstruct laboratory-measured phy- 
toplankton in vivo absorption spectra has been considerably 
improved, with the results of the present work demonstrating 
coefficients of determination across the visible spectrum at 
greater than 0.9 for several of the inversion techniques em- 
ployed. The solutions obtained by reconfiguring the equations 
to invert across the spectrum of a single sample to obtain 
pigment estimates are less successful. Of the 12 inverse model 
solutions carried out to estimate pigment concentrations from 
phytoplankton absorption spectra, none had all good (r 2 > 0.6) 
retrievals, for either the U.S. east coast or global ocean data 
sets. The best performer using the more regional U.S. east coast 
observations is from using the SVD-NNLS method, which 
uses the SVD algorithm to obtain pigment-specific absorption 
spectra and then the NNLS algorithm to solve for the pigments 
using the observed phytoplankton absorption spectra. The 
results are significantly improved for most of the inversions if 
the absorption spectra are first ‘normalized’ to the expected 
675 nm absorption levels. 

[ 45 ] Most pigments show either consistently good (chlo- 
rophyll a, chlorophyll c, peridinin, fucoxanthin, diadinox- 
anthin, and carotenoids) or poor (chlorophyllide, both 
fucoxanthin, diatoxanthin, zeaxanthin, and lutein) retrievals. 
There are a number of factors that are likely contributors to 
this. The first is related to the ESD approach that is used to 
generate pigment-specific absorption. A number of the 
pigments are correlated (Figure 4), adding complexity to the 
inversion solutions. In addition, the solutions explicitly 
contain the influence of the pigment package effect. Both of 
these add more degrees of freedom to an already overly 
complex set of observations. Also, the spectral inversions 
are done at a spectral resolution of 1 nm, which is of the 
order of spectral distance between several absorption spectra 
(Figure 1). Results from twin experiments to test the 
development of these inversion techniques (not shown) 
show that inverse model solutions are sensitive to the 
spectral resolution used and the level of errors in the mea- 
sured phytoplankton absorption spectra. 

5. Conclusions 

[46] Phytoplankton absorption spectra and pigment 
observations can be used with an inverse modeling technique 
to support a second inversion method that estimates phyto- 
plankton pigment concentrations directly from phytoplank- 
ton absorption spectra. This capability can be directly inserted 
in present hyperspectral remote sensing inverse models to 
expand the pigment products being recovered. The lack of 
measurements of MAA’s, phycobiliproteins and other pig- 
ments impacts the capability of this technique to retrieve 
pigment-specific absorption spectra and estimates of pigment 
concentrations. In addition, the pigment estimations are less 
accurate for those pigments whose combined absorption le- 
vels and concentrations have a low influence on the final total 
absorption spectra. Finally, the most accurate model for 
estimation of phytoplankton total absorption, which possibly 
accounts for pigment package effect and a missing absorption 
term, is less able to be inverted to estimate phytoplankton 
pigments that the more robust and straightforward ESD 
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method that calculated pigment-specific absorption ‘modes.’ 
The SVD inversion is the most robust for inversion applica- 
tions, even with its problem in not being able to guarantee 
nonnegative absorption or pigment values. The inversions 
show good retrievals for a number of phytoplankton pigments 
from water samples obtained within a region of the U.S. East 
Coast. The capability to obtain pigment retrievals is lower 
when using a larger global ocean data set. It is likely that this 
method can only be best applied in a regional setting. 
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